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I.  INTRODUCTION 

The  purpose  of  this  report  is  to  present  a number  of  re- 
sults regarding  errors  in  radar  measurements  produced  by  propa- 
gation through  the  ionosphere.  The  emphasis  is  on  the  dependence 
of  these  errors  on  radar  frequency,  and  on  the  possibility  of 
determining  and  removing  the  errors  by  making  simultaneous 
measurements  at  two  radar  frequencies. 

The  application  which  motivates  this  collection  of  results 
is  primarily  to  data  taken  by  the  ALTAIR  radar,  which  operates 
simultaneously  at  VHF  and  UHF  frequencies.  For  these  frequen- 
cies, and  particularly  at  VHF,  errors  in  metric  data  induced  by 
the  ionosphere  are  often  large  in  comparison  with  errors  from  all 
other  sources.  Thus  the  degradation  of  metric  data  caused  by 
propagation  through  the  ionosphere  is  significant,  and  the 
measurement  of  this  error  is  highly  desirable. 

The  use  of  two-frequency  measurements  to  determine  iono- 
spheric corrections  is  by  now  a well-established  technique  , ^ ^ 
but  until  recently  these  methods  have  not  been  applied  routinely 
to  data  from  the  KREMS  radars.  Efforts  toward  correcting  the 

metric  data  have  focused  primarily  on  computing  the  correction 

(2) 

using  model  ionospheres  and  ray-tracing  calculations.  It  is 
the  purpose  of  this  report  to  show  the  feasibility  of  determining 
these  corrections  by  direct  measurement  using  current  ALTAIR 
capabilities.  The  techniques  are  currently  being  used  in  post- 
mission  analyses  with  gratifying  success. 

The  ionospheric  effects  discussed  here  are  those  on  radar 
range,  elevation  angle,  coherent  doppler  velocity,  pulse  shape, 
and  wave  notarization.  With  the  exception  of  the  elevation 
correction,  these  results  are  well  known.  They  are  presented  for 
convenience,  ana  in  order  to  cast  them  info  form  which  demon- 
strates their  applicability  tu  ALTAIR  data  fo  lections  of  the 
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azimuthal  angle,  effects  due  to  ionospheric  disturbances  , and 
tropospheric  refraction  are  not  considered  here. 

2 . RANGE  ERROR 

The  apparent  range  R of  a body  from  the  radar,  as  measured 
by  time  delay,  is: 

R = c / ds/v 
Jp  g 

where  v is  the  group  velocity  and  the  subscript  p indicates  that 
g 

the  integral  is  one-way  over  the  ray  path. 

For  ionospheric  propagation,  v^  = cn , where  n is  the  re- 
fractive index.  If  the  effect  of  the  earth's  magnetic  field  and 
electron  collisions  can  be  ignored,  n is  given  by 

n2  = 1 -2X 

with 

X = 40 . 2N/ f 2 

-3 

where  N is  electron  density  in m and  f is  radar  frequency  in 

Hz.  For  VHF  frequencies  or  higher,  X « 1 and  terms  of  order 
2 

X can  be  ignored.  Thus 

R = i » Id sd  + X) 

•'p  n -'p 

The  effect  of  path  curvature  on  apparent  range  can  be  shown  to 
be  0(X2)  so  that  / ds  » R . , the  true  range,  and 

t 

LR  = R - R = / dsX(s)  (1) 

t J p 

is  the  range  error.  Rewrite  R as 

R = Rt  + B/f2  (2a) 

whe  re 

B = 40.2  /dsN(s)  (2b) 

If  R is  measured  at  two  different  radar  frequencies  ^ and 
f the  measured  ranges  differ  by 
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ar 
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R. 


) 


also 


-R!  R1  ~ Rt  Bfl 


_? 


Eliminating  B: 


R1  - Rt  = g • (R2  - R1)  <3> 

where 

g = [il/i2)2  - l]'1 

Thus  the  range  difference  AR12  can  be  used  to  find  the  range 
error  AR^,  the  true  range  R^ , and  the  apparent  range  at  any  other 
frequency  f„.  In  particular,  for  ALTAIR  we  can  set  f1  = 415  x 10  , 
f = 155,5  x 10  , giving  g = 0.16333.  Using  this  value  in  (3) 
with  the  UHF  - VHF  range  difference  gives  the  UHF  range  error 
and  also  VHF  error. 

3.  PULSE  BROADENING 

Once  the  range  error  is  computed  by  two-frequency  measure- 
ment, the  amount  of  spreading  of  either  pulse  due  to  ionospheric 
dispersion  can  be  found,  provided  that  the  undispersed  pulse 
shape  is  known.  If  S(t)  is  a chirped  pulse,  and 

S(oj)  = jfx  dt  S ( t ) e-la)t 


is  its  fourier  transform,  then  the  signal  received 
can  be  written  as 

CC  -W  ^ 

A(t)  = / — e140t  S (oo)  P(oo)  S*  (oo)  W (W) 

2rr 


where 

P(oo) 


by 


the  radar 
(5) 
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is  the  transfer  function  for  propagation  through  the  ionosphere. 
Here  S*(oo),  the  complex  conjugate  of  S,  is  the  matched  filter, 
and  w(w>  is  any  additional  weighting  applied  to  reduce  sidelobes. 
Expression  (5)  can  also  be  written  as 

A(t)  - f Ao(a))PC^ 

“'-CD 

where  _ ^ 

A (go)  = S(W)  2 W(w)  is  the  transform  of  the  compressed 
pulse.  Expanding  the  refractive  index  n as  before: 

A(t)  _ f dw  eiw(t  - ~)  A ((0>  e1  Z£  (6) 


with  E = 8tt  B. 


The  spectrum  A^Coo)  is  a peak  located  at  the  radar  angular 
frequency  u>  , so  that  Aq(go)  = A^u)  with  u=w-W0,  and 

a,  . r 2 E | (7 

iou  t f ciu  A ( * exfut-  + u “3~J 
A(t)  = e o + J_m  A^u)  e 


with 


2R 

t+  = t - — + 


Here  we  have  expanded  of1  in  powers  of  u and  retained  0(u  ). 

If  the  compressed  undispersed  pulse  AQ(t)  is  approximated 

by  a Gaussian  : 2 


A (t)  oc  e ia)ot  e"l^ 
o 


then 


A^(u)  oc  e 


-(tu) 


Inserting  this  form  in  (7),  we  find  2 

A ( t ) « e1W°t+  (a2-  iO”2  e 4(a*  - ie) 


with 


a = oo  ' and  •-  = E/(oo  c). 
wo  ° 
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The  pulse  shape  in  terras  of  radar  cross  section  is  pro- 
portional to 


|A(t)| 


2 , 4 2s-k  - 

“ (a  + e ) e 


, + ,2  2 
(%t_)  a 


“5 2 

a + € 


(10) 


From  this  it  follows  that  the  effect  of  ionospheric  dis- 
persion is  to  increase  the  width  and  decrease  the  peak  height 
by  a factor 

F = [1  + Y2  (11) 

2 2 

with  Y = e /a  = L R.  For  a given  frequency  go  and  com- 

U)0CT  ° 

pressed  pulse  width  t,  the  factor  F(aR)  determines  the  shape  of 

the  dispersed  pulse.  For  ALTAIR  VHF,  fQ  = 155.5  MHz  and  a 3db 
width  of  0.195u,sec  results  in  t 0.083^scc  and  Yy  = 0.986£Rv, 
where  ARv  is  VHF  range  error  in  kilometers.  For  UHF,  fQ  = 415 
MHz  and  a 3db  width  of  O.luis  implies  t « 0.042u-s  and  yu  = 1.42<1Ru 
with  AR  the  UHF  range  error  in  kilometers.  Thus,  a measured 
AR  range  difference  of  2 km  corresponds  to  a spreading  factor 
F - 2.5  at  VHF,  and  F = 1.05  at  UHF.  Since  2 km  is  a fairly 
typical  range  difference  observed  in  tracking  bodies  above  the 
peak  of  the  ionosphere,  it  can  be  concluded  that  at  these  target 
heights  the  VHF  pulse  will  often  be  broadened  by  a factor  of  two 
or  more,  while  the  dispersive  effect  on  the  UHF  pulse  is  negli- 
gible . 

The  question  remains  as  to  whether  the  Gaussian  approxi- 
mation of  the  compressed  pulse  is  a good  one.  Figure  1 shows 
a measured  VHF  long  chirp  compressed  pulse  shape,  together  with 
a fitted  Gaussian.  The  pulse  shape  data  are  from  the  test  ALTO 
319B1  of  31  March  1971.  The  data  points  were  generated  by  per- 
forming a.  .range  scan  on  a calibration  sphere.  The  gaussian  was 
fit  by  eye  to  the  data  points,  which  determined  the  value  of 
given  above.  The  measured  pulse  is  seen  to  be  slightly  skewed 
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(0P) 


at  the  top  of  the  leading  edge,  but  over  most  of  the  main  lobe  the 
agreement  is  very  good. 

If  greater  accuracy  than  that  provided  by  the  Gaussian  fit 
is  required,  pulse  shape  dispersion  can  be  computed  without 
approximation.  This  is  done  by  calculating  numerically  the 
convolution  integral  of  the  measured  pulse  shape  with  the  time 
transform  of  the  propagation  factor.  However,  the  reasonably 
good  fit  by  the  Gaussian  spectrum  indicates  that  the  factor  F 
given  by  (11)  provides  a sufficiently  accurate  measure  of  the 
pulse  broadening. 

The  pulse  spreading  has  significance  for  both  radar  cross 
section  (RCS)  and  range  measurements.  If  the  VHF  RCS  of  a hard 
body  is  measured  as  the  return  from  the  peak  of  the  pulse,  then 
the  measured  RCS  will  be  too  low  by  a factor  of  F.  With  RCS  in 
dbsm , the  error  is 

'■'°db  - -10  loB10  F (l2) 
From  the  numerical  estimates  above,  it  follows  that  pulse  spread- 
ing can  cause  errors  in  measured  RCS  as  large  as  5db  at  VHF. 

Similarly  the  range  resolution  of  the  VHF  pulse  is  taken 
to  be  the  3db  width  of  29m,  as  in  Figure  1.  The  effect  of 
propagation  through  the  ionosphere  is  to  multiply  this  number  by 
F and  thus  degrade  the  resolution  by  this  factor.  The  effect 
should  be  taken  into  account  when  measurements  depending  on  the 
resolution  length  are  analyzed. 

One  example  of  this  dependence  is  the  measured  RCS  of  chaff 
clouds.  Although  pulse  spreading  causes  the  RCS  of  an  individual 
dipole  to  decrease  according  to  (12),  the  range  resolution  cell 
increases  in  length  by  a factor  of  F,  so  that  on  the  average,  F 
times  as  many  dipoles  contribute  to  the  return.  Thus  except  for 
very  sparse  chaff  clouds  (less  than  one  dipole  per  resolution  cell) 
and  except  for  positions  near  the  edge  of  the  cloud,  the  average 
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chaff  RCS  should  be  unaffected  by  the  pulse  broadening.  Since 
the  RCS  of  a hard  body  immersed  in  chaff  is  decreased  by  (12) , 
the  signa 1-to-noise  ratio  (S/N)  of  body  to  chaff  RCS  is  also  de- 
creased by  an  amount  (12).  Therefore  the  pulse  spreading  can 
cause  losses  of  5db  in  S/N  when  searching  at  VHF  for  a hard  body 
in  a chaff  cloud. 


4.  ELEVATION  ERROR 

A.  Flat  Ionosphere  Approximation 

Once  AR  is  known  for  a given  radar  frequency  f > 100  MHz , 
the  error  in  elevation  angle  can  be  computed  without  ray-tracing. 
To  see  this,  consider  first  a flat  ionosphere  with  electron 
density  N(h)  a function  of  height  only. 

Referring  to  Figure  2 , Eq  is  the  apparent  elevation  , t 
the  true  elevation,  and  6 = Eq  - t is  the  elevation  error.  The 
ground  range  xt  can  be  represented  in  two  ways: 

x = r cos  t = / dx  = / ds  cos  E(h)  (13) 

t t J P JP 

Snells'  Law  is 


n cos  E = const. 

= cos  E 

o 

since  N(h)  is  assumed  to  vanish  at  the  radar  (N(0)  = 0)  . 
bining  these  two  expressions: 


R 


t 


COS  T 

cos  E 

o 


r ds/n 
JP 

R (the  apparent  range) 


Com- 


If  the  error  6 is  small,  we  can  write 

cos  t = cos(E  - 6)  « cos  E + 6 sin  E 

O O U 

so  that 

1 + 6 tan  Eo  - 
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where  we  have  replaced  by  R in  the  last  expression,  intro- 
ducing an  error  of  order  . Equation  (14)  allows  one  to  find 
the  true  elevation,  given  only  the  measured  range  and  elevation, 
and  the  range  error.  It  is  valid  to  the  extent  that  the  earth 
can  be  thought  of  as  flat,  and  assuming  6 is  small  . In  particu- 
lar, the  approximation  of  cos  t requires  6<<:1  and  6<<2  tan  Eq. 

For  radar  frequencies  f > 150MHz  and  Eq  > 1 , these  inequalities 
are  virtually  always  satisfied. 

Equation  (14)  was  checked  against  the  results  of  ray-tracing 

computation  of  6 . A re-entry  trajectory  and  density  profile 
N(h)  were  chosen  arbitrarily,  and  the  ray-tracing  program  com- 
puted the  corrections  ^R  and  6 for  f = 155.5  MHz.  The  com- 
puted 6's  agreed  with  (14)  to  within  about  30%  during  rater 
portions  of  the  trajectory  (R  < 1000km)  , within  about  50%  at 
2000km  range  and  within  a factor  of  2 at  2750km.  It  is  therefore 
desirable  to  recalculate  (14)  for  the  spherical  case. 

B.  Spherical  Case 

Referring  to  Figure  3,  the  differential  elements  are  re- 
lated by 

rd9  = os  cos  E(h) 
dr  = ds  sin  E(h) 

Snell's  Law  for  the  spherical  case  is 

nr  cos  E = r cos  E 
o o 

Using  these  relationships  we  can  again  construct  an  identity 
from  which  6 can  be  computed.  Instead  of  ground  range,  it  is 
more  convenient  in  this  case  to  use  the  geocentric  angle  0 : 


where 

as 


re 

r . cos  F. 

6 = J 

de  = 
o 

r 

= j ds  = 

Jp  r 

/ \2 

1 

f dr 

( nr  ) “■> 

= J 

r r 

A 

\rocos  Eo/ 

r is 

the  earth  radius.  T1 

r 

L 


dr 


cot  E 


(15a) 


1 

~2 


o 


6 = tan 


-1 


R , + r sin  t 
t o 


r cos 
o 


(15b) 


We  then  proceed  as  in  the  flat  case,  equating  these  two 
expressions  for  e,  expanding  and  keeping  terms  of  first  order 
in  X and  in  6-  The  result  is 

6 = cot  E (l  + n) 


du 


X (u) 


o 


(1  + Bu) 


A 

R = R / duX(u) 


(16a) 

03b) 


where 


a nd 


R 

” r sin  E 
o o 


u = S/R 

Here  S is  the  length  along  the  straight-line  path  to  the  body 
being  tracked.  Equation  (16b)  is  equivalent  to  (1),  since  the 
path  difference  introduces  an  error  of  order  X . 

C.  The  Sharp  Ionosphere  Approximation 

In  order  to  again  express  6 in  terms  of  AR,  we  assume  as  a 
first  approximation  that  the  ionosphere  is  sharply  peaked  at  a 
height  H.  corresponding  to  a value  u = up.  Then  X(u)  is 
essentially  a 6-function,  and  we  have 
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6 


du  X ( u) 


cot  E 


1 + B 
( '■  + Bu  )' 


/„ 


i + e 


(1  + Pu  ) 
p 


4-1?  R 

cot  E — 
° R 


w he  re 


(1  + eu  )2  - 1 + H(gr°  ; H) 

P r “sin  E 

o o 


(17) 


(18) 


Equation  (18)  defines  u 

P 

height  H,  i.e.  for  u^  < 1.  For 

u = 1 in  (17) . 

P 

Equation  (17)  reduces  to 
flat-earth  result  (14)  is  valid 
In  the  opposite  limit,  R 
cos  E„ 


W 


r.R 


as  long  as  the  body  is  above  the 
lower  altitudes  we  can  set 

(14)  as  B - o,  showing  that  the 

for  R « r sin  E . 

o o 

» r sin  E , (17)  reduces  to 

o o ’ 


This  approximation  holds  early  in  an  observed  re-entry, 

when  E is  small  and  R is  large.  In  this  case,  cos  E 1,  and 
o 

if  we  arbitrarily  choose  an  effective  height  H = 400km  for  the 
ionosphere,  we  arrive  at 


1.25 


uR. 


mr  ~ km 

where  6 is  in  milliradians  and  R is  in  kilometers.  This  ex- 
plains tie  rough  equivalence  belween  6 and  LR  which  has  been 
observed  in  comparing  ALTAIR  and  ALCOR  metric  data  from  simul- 
taneous tracks  on  the  same  target.  The  fact  that  the  equivalence 
does  not  always  hold  is  also  explained,  since  it  is  only  valid 

for  R » r sin  E . 

o o 

Equation  (17)  was  tested  against  the  ray-tracing  program 
using  a re-entry  trajectory,  a satellite  trajectory,  and  several 
ionospheric  profiles  in  various  arbitrary  combinations.  It  was 
found  that  with  a proper  choice  of  H,  6 agreed  with  the  ray- 
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tracing  result  to  w .thin  2 or  3%  for  body  altitudes  sufficiently 
highei*  than  H (about  100km  higher).  At  lower  altitudes,  the 
agreement  typically  worsened  to  10  to  15%.  To  improve  this 
situation  it  is  necessary  to  include  the  effect  of  non-zero  width 
of  the  ionospheric  peak. 


I).  The  Gaussian  Ionosphere  Approximation 

A rough  model  of  the  peak  is  a Gaussian: 
[h  - H 12 


w 


N(h)  = N e 
o 


hich  translates  near  the  peak  to 


X(u)  = Xqc 


u - u 


with 


a = oh/  Jr  Sin  Eq  (1  + 0upCos2Eo)| 

Using  this  X(u)  in  (16a),  expanding  and  keeping  terms 
first  order  in  o: 


^ T,  £R  1 + P 
6 = cot  E — rf 

° R (1  + 0u  )2 
P 


1 + 


Bo 


-A 


1 + Bu 


p L 


e y dy 


1 — u 


where  A = 


and 


fk  e-^dy  * J A"  |l  + =rf  A| 


(19) 

of 


(20) 


In  equation  (20),  up  is  given  by  (18)  even  for  up  > 1.  The 
effective  width  ah  of  the  ionosphere  varies,  but  after  examining 
several  profiles,  a value  ah  = 150km  was  chosen  as  a reasonable 
average  and  was  used  for  all  cases.  With  fixed,  (20)  is  a 
one-parameter  formula  for  the  elevation  error. 

Equation  (20)  was  also  tested  against  the  ray-tracing 
program.  To  do  so,  a re-entry  trajec  ory  was  generated  by  the 
NRTPOD  (Non  Real  Time  Precision  Orbit  Determination)  program,  using 
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a state  vector  from  a recent  KREMS  test.  This  trajectory  was 
used  as  input  to  the  ray-tracing  program,  which  computed  by 
iteration  an  apparent  trajectory  for  a radar  frequency  of  155.5 
MHz,  using  one  of  three  different  ionospheric  profiles  taken  at 
Jicamarca,  Peru.  The  apparent  range  and  elevation  and  the  range 
error  were  then  used  in  (20)  to  find  6,  which  was  compared  with 
the  elevation  error  from  the  ray-tracing  program.  The  same  pro- 
cedure was  followed  with  a satellite  trajectory  using  a pass  with 
a maximum  elevation  of  81°  and  minimum  range  of  800km.  The 
effective  height  H of  the  ionosphere  was  varied  in  (20)  to  find 
the  best  fit  for  each  profile  and  also  the  sensitivity  of  the 
fit  to  errors  in  H.  The  best  value  of  H appeared  to  be  either 
the  peak  height  of  each  profile  or  slightly  higher,  depending  on 
the  asymmetry  of  the  peak.  An  error  AH  in  estimating  H produces 
an  error  46  in  6,  where  roughly 

A6  ^ (21) 

6 H + \ v sin2E 

£ o o 

Thus  the  percentage  error  in  6 is  less  than  that  of  H. 

Tables  I - VIII  show  the  results  for  two  of  the  three 

* 

ionospheric  profiles,  shown  in  Figure  4.  and  both  trajectories. 
Delta  EL  is  the  elevation  error  computed  by  ray-tracing,  and  Delta 
EL2  is  6 from  equation  (20).  The  relative  error  is 

AEL  - 6 
AEL 

It  is  seen  from  able  I-IV  that  by  choosing  H carefully, 
one  can  predict  the  el  it  ion  error  to  within  3^  over  the  whole 
trajectory.  This  can  be  taken  as  verification  of  equation  (20). 


* 

The  names  NORM40  and  LOW2-31  simply  label  the  profiles  and  have 
no  significance. 
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TABLE  I 


noowuo 


TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 'i  50  km 


nKVSTTY  PP^NTRY  TM.TPCTnpv  1c6.6f»  KH’’ 


PA  MOP 

OP  I T A 
R A NOP 

PLPV  A’T’TOR 

1674.14 

1.7007 

0.2210 

152  0. 74 

1 .7046 

1. 1616 

1164. °8 

1 .6010 

2.6011 

1707.00 

1,6  706 

1.  6160 

1046.75 

1.6620 

4.7654 

7304. 1 0 

1 .6201 

6.8011 

77  10.  ? 1 

1.6011 

7.0101 

2561 . 05 

1.6574 

0. 1416 

2781.05 

1.6177 

«.  2^27 

270°. 56 

1. 4714 

10. 1070 

2034. 57 

1.4220 

11. 4000 

1 056. 04 

1.1630 

12.6110 

1676.4? 

1.2011 

13. 7210 

1401.21 

1.1040 

14.0221 

1107. 10 

1.0472 

15.0050 

1110.20 

0,.  0 767 

16. 0500 

026.54 

0,6060 

17. 0*20 

712.26 

0.1071 

10.0427 

516.76 

0.0161 

10.6004 

140. 00 

0.00  10 

10. 1670 

DEI.T  A 

TELTA 

PIT  A'T’T V ^ 

PL. 

PT.l 

FpnnP 

0,  14046 

0,11710 

0.07110 

0.  14426 

0.  140O0 

0.o?ni 

0 . X 4 660 

0.  14147 

0.02116 

1.  14750 

0.  14405 

0.0 1046 

0.  14714 

0.  14516 

0.01406 

0.  146  11 

0.  1445? 

0.01000 

0.  14414 

0.  V411  3 

0.00606 

0.14167 

0.14110 

0.OO141 

0. 11006 

0.  11080 

0. O0O60 

0.  13610 

0.11611 

0.00107 

0.  11145 

0.1116? 

0.O01OR 

0„ 11000 

0. 11001 

0.00  0 10 

0.  12004 

0. 12707 

0. O01 10 

0. 12440 

0.12416 

0.00  100 

0.  11741 

0. 11714 

0. 00070 

0.  10104 

0.  10110 

0 ,002c1 

0 07304 

0. 07101 

0.001 10 

0.03410 

0.01100 

0 . O 1 40  1 

0.00000 

0.00860 

0. 01465 

0. 00040 

0.  00040 

0.006OQ 
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TABLE  II 


KinpMon 


TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 350  km 


DENSITY  SATET  TPA.IECTOBY  155.50  "87 


DETTA  HFLTA 

PANGF  PANGE  ELEVATION  EL. 


3306. 21 

1.6079 

311R.47 

1.6907 

2930.37 

1.6677 

2743.  61 

1.6280 

2656.96 

1.571R 

2371.25 

1.5003 

21R6.07 

1.4156 

2004. 63 

1.3203 

1R24.R6 

1.2175 

1648.99 

1.1101 

147R.U6 

1.0013 

1315.53 

0.R943 

1163. 5R 

0.7927 

1027.66 

0.7009 

Q 1 5.  07 

0.6245 

035.44 

0.5704 

79  R . 75 

0.5455 

RIO. 00 

0.5539 

069. R3 

0.5943 

966. OR 

0.6607 

1092.14 

0.7459 

1236.73 

0.0437 

1394.62 

0.9490 

1561.67 

1.0578 

1735. 10 

1.1668 

1 °1  3 . 05 

1.272R 

2094.22 

1.3724 

2277.7  1 

1.462R 

2462. R5 

1.5409 

2649.13 

1.6045 

2036.  10 

1.6517 

3023.73 

1 ,6R  19 

3211.52 

1.6957 

-0.0839 

0. 13875 

1,6412 

0. 14539 

3. 4753 

0.  14873 

5.4406 

0. 148  24 

7 *>663 

0.  14309 

9.0905 

0.  13611 

12.4637 

0.  12561 

15 . 3533 

0.11326 

1R.6496 

0.  09990 

22.4765 

0. 08620 

27.0039 

0.07267 

32.4654 

0.05965 

39. 1729 

0.  0472R 

47.5058 

0.03560 

57.RC63 

0.0245P 

69.9886 

0.01427 

81.  1651 

0.00609 

76.2040 

0. 00962 

63, 5708 

0.01936 

52. 3434 

0.03010 

43.0395 

0. 04152 

35.5547 

0.05363 

29. 5036 

0. 06644 

24.5  351 

0. 07985 

20. 3767 

0.09362 

16.B2R6 

0.10731 

13.7453 

0.  120  30 

11.0207 

0.  131B1 

R. 576° 

0.  14094 

6. 3554 

0.  146R9 

4. 3124 

0.  14906 

2. 4142 

0.  14729 

0.6356 

0.  14190 

DEL^A 
EL  2 

0. 13488 
0.14151 
0.  14530 
0. 14596 
0.  14294 
0.13546 
0.  12699 
0, 115 >6 
0. 10206 
0.0RR22 
0.07435 
0.06091 
0.  04 P 1 4 
0.0:1613 
0.0/4RR 
0. 01441 
0.00614 
0.00071 
0..  01  057 
0.03051 
0.04220 
0. 05460 
0.06702 
0.  OR  172 
0.09574 
0. 10941 
0. 12100 
0. 13265 
0. 14056 
0.  14514 
0.14609 
0.14353 
0. 1379  3 


BET  A^IVF 
EPPOP 

0.02703 

0.02670 

0.02243 

0.01540 

0.00661 

0.00250 

0.01090 

0.01751 

0.02166 

0.02343 

0.02306 

0.0210O 

0 . 0 1 R 25 

0. 01400 
0.01210 
0.00999 
0.00900 
0.00030 
0.01007 
0.01362 
0.01661 
0.01973 
0.02220 
0.02339 
0.02265 
0.0195R 
0.01406 
0.00636 
0.0026R 
0.01191 
0.01903 

0.02652 

0.02001 
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ABLE  III 


TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 385  km 


Tnw?-?1  OPNSTTY  OFF  NTR Y T R A.1  DRY  1^5.50  WM? 


nPT.TH 

DPT.  7*  A 

OFT  TA 

PPT  AT  T V p 

ramoe 

RUMOR 

PLPV A^TOM 

PL. 

PL? 

pnnOR 

1873.05 

0.4174 

0.1115 

0. 030  18 

0.03075 

0.00??? 

3819.45 

0.4180 

1.2503 

0.03101 

0.03111 

0. no? 07 

3383.70 

0.4111 

2. 3381 

0.01181 

0.03177 

0 . 00  ?83 

3205. 73 

0.4007 

3.5103 

0.01104 

0.03200 

0.00438 

3045.50 

0.4020 

4.8501 

0.01207 

0.03772 

0. O0878 

23R2.05 

0.1058 

5.7730 

0. 03100 

0.03214 

0.00753 

2710.03 

0. 1370 

8.O053 

0.03181 

0.01137 

0. 003?? 

2550.87 

0.1783 

3.0310 

0.03113 

0.01144 

0.0031? 

2330. fll 

0.1850 

0.1544 

0.03085 

0.03035 

0.OO885 

2200.44 

0.1512 

10.2753 

0.03004 

0.03014 

0. 00?40 

2031.45 

0.1140 

11. 3040 

0.02014 

0.02010 

0.00150 

1055.70 

0.1150 

12.5107 

0.02352 

0.  02310 

0. 00770 

1575. 42 

0.2003 

13.8211 

0 . 027  4 0 

0.07711 

0. 01383 

1402.29 

0.2504 

14.7733 

0. 07803 

0.07588 

0.018  14 

1->08. 38 

0.2137 

15.3174 

0.02338 

0.02354 

0.01344 

1117.82 

0.1828 

18. 3753 

0.01033 

0.0137? 

0.003  I? 

°28 .13 

0.004 1 

17.3028 

o. 01145 

0.01740 

0 .00304 

7?2. 10 

0.0271 

13.3111 

0.00435 

0.00473 

0.01318 

*38.21 

0.0030 

10. K017 

0.0007? 

0.00071 

0 . 01007 

•*40.30 

0.0003 

10.  3887 

0.00003 

0.00010 

0.27307 
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TABLE  IV 


TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 385  km 

T o v 2 - 1 1 OF V 6 T f Y PA^FT  ^ F A.lFC'T,nPY  1*6. *0  H 


0 0 T o’  A 

0 FLT  A 

DELTA 

8PT  A^T VP 

PA  VOP 

p AMOK 

FLFV  ATTON 

FT. . 

FT,  7 

FP8H8 

7704.8? 

0.4036 

-0.1833 

0. 02043 

0.07876 

0.81681 

7 117.10 

0.4074 

1. 6267 

0.  03007 

0.07047 

0.8147. 

2020.  f-n 

0. 3070 

3.3684 

0.03184 

0.07148 

0.01107 

?742. 37 

0. 380p 

*. 3744 

0. 03204 

0.07188 

0.8048? 

?666.76 

0.3781 

7. 4830 

0.03148 

0.03158 

0.0071? 

2 ’7n.  1 1 

0. 3620 

Q.7846 

0. 03070 

0. 03056 

0.01165 

7186.86 

0.  3444 

17. 3664 

0. 02830 

0. 07886 

0.0188q 

2007 .*7 

0.3231 

16. 7668 

0. 07681 

0. 02668 

0. 8 2 6 8 

1P23.0/; 

0. 2°07 

18.6720 

0 . 02370 

0.  02788 

0.07888 

ifup.  16 

0. 7748 

72.  4106 

0.07031 

0.02"76 

0.07167 

1077.71 

0.2402 

26.0486 

0.017  35 

0. 01^80 

0 . 8708  7 

1314.  86 

0.7736 

37.4202 

0.0144? 

0. 0 148  3 

0.07P61 

116?. oo 

0.  10-30 

30. 1371 

0. 01 166 

0.01183 

0.0761? 

1027.13 

0.1766 

47. 4780 

0. 00877 

0.00886 

0.07088 

Q 1 4 . f 0 

0.1677 

67.7878 

0 . 00600 

0. 00670 

0.01778 

876.01 

0.  1447 

60.0670 

0.00356 

0.00360 

0.01487 

708.^4 

0.1300 

8 1.  160  8 

0. 00167 

0.00164 

0.01348 

R1 0.40 

0.  1401 

76. 1 c6  8 

0.  00240 

0.00243 

0.0148? 

860.30 

0.  1601 

63,  6867 

0.00481 

0 . 004  P 0 

0.01586 

066.40 

0.1666 

82.  3707 

0.00744 

0.00758 

0. 01880 

1001.88 

0.1878 

43. 0087 

0.0  10  18 

0.01041 

0 . 0??8R 

1236.10 

0.2113 

38.8141 

0.01707 

0. 01337 

0. 07668 

1 303.01 

0.7366 

20.4532 

0.01508 

0.01643 

0.07867 

1660.87 

0.7624 

24.4741 

0.01803 

0.  01867 

0.03176 

1734.22 

0.7878 

20. 3080 

0. 07188 

0.02766 

0.03061 

1Q12.00 

0.3171 

16. 7460 

0. 07471 

0.02538 

0.07758 

2003. 10 

0.3346 

13.6*37 

0. 02777 

0. 02788 

0. 0??10 

2276.60 

0.3843 

10.0183 

0.  02043 

0.07886 

0.01464 

246 1. 68 

0. 3711 

8 . 4670 

0.03101 

0. 03120 

0. 00612 

2647. 01 

0. 3846 

6. 2404 

0.03180 

0.0318? 

0.00733 

2834.04 

0 . 3 04  1 

4.  1083 

0.07201 

0.01171 

0.00866 

3022. 48 

C.  40 ')1 

7.7803 

0. 071 78 

0.0308  7 

0.0146? 

3 210,,  37 

0.4076 

0. *730 

0.03010 

0.07860 

0.0167? 
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Tables  V-VIII  show  that  if  c ne  makes  an  error  of  50km  in  esti- 
mating H,  then  (20)  gives  a result  correct  to  within  15%  over 
the  entire  trajectory.  The  maximum  absolute  error  of  6 in  this 
case  is  about  0.02°  at  VHF , compared  to  a total  elevation  error 
of  0.14°.  Thus  most  of  the  elevation  er’*or  due  to  ionospheric 
refraction  can  be  removed  by  use  of  (20)  , even  if  the  height  H 
is  not  well  known. 

E.  Estimation  of  the  Ionospheric  Height  H 

With  the  width  parameter  a fixed  in  equation  (20) , four 
inputs  are  required  by  this  expression  in  order  to  compute  ele- 
vation error.  These  inputs  are  the  measured  range,  elevation, 
and  two-frequency  range  difference,  and  the  height  H.  The  method 
of  choosing  H can  be  different  for  real-time  vs.  post-mission 
corrections,  and  may  differ  for  different  types  of  missions. 
Possible  methods  include  the  following: 

1.  Set  H - 375km  (or  some  other  constant  average  height).  The 
true  H is  unlikely  to  be  less  than  300km  or  greater  than  450km, 
so  the  maximum  error  is  roughly  75km,  corresponding  to  an  error 
in  6 of  less  than  25%.  This  is  the  simplest  solution  and 
probably  adequate  for  real-time  estimation. 

2.  Make  use  of  a priori  knowledge  of  diurnal  and  seasonal 
variations  in  the  height  of  the  ionosphere.  Use  as  input  a value 
consistent  with  time  o_  day,  time  of  year,  location  in  the 
sunspot  cycle,  etc.  Tables  or  curves  of  the  variation  of  H could 
be  generated. 

3.  Use  measurements  of  the  ionosphere.  If  electron  density 
profiles  taken  just  before  a mission  could  be  obtained  from  an 
ionosonde  located  at  KREMS  or  at  some  point  near  the  ground  trace 
of  the  trajectory,  the  effective  height  could  be  estimated. 

4.  Differentiate  the  range  error  "R.  If  the  look  angle  (the 
angle  between  the  target  velocity  vector  and  radar  line-of-sight) 
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TABLE  V 

TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 475  km  Instead  af  385  km 


LOW  2-3  1 DENSITY  REENTRY  TRAJECTORY  1^5.50  MU Z 


FANCE 

DELTA 

RANGE 

FIE VATICN 

3 6 73  . C 5 

C.4174 

C.1115 

3 519.  46 

C. 4 160 

1.2503 

3363.70 

0.4131 

2. 3861 

3:0  5.'i  3 

C.40e7 

3.5153 

3045.50 

C . 40  29 

4.C5C1 

2662.95 

C . 3956 

5.7789 

271 G. 03 

C.3070 

6.9C58 

255C.67 

C.3760 

8.C310 

233C.63 

C.365C 

9.  1544 

2 2C  6 . 44 

C.3512 

1C. 2758 

2033.45 

C.3349 

11. 3949 

1655.79 

0.3150 

12.5107 

1675.42 

C .2903 

13.6213 

1492.29 

C . 2 5 Q4 

14.7238 

1306.36 

C .2187 

15.8124 

1117.62 

C.  1626 

16.8758 

5 2 6.  13 

0.C941 

17.8526 

732.10 

C.C271 

18.8121 

536.21 

C.OC30 

19.5012 

34 C . 60 

C.OOC3 

19.3667 

DELTA 

DELTA 

REI ATI VE 

EL. 

EL2 

EFECF 

0.03015 

0. C2667 

0.11540 

C „ 0 3 1 0 1 

0.02744 

0.11519 

0.03161 

0.028C1 

0.11361 

0.03194 

0. 02839 

11092 

C.C3202 

0. C 2850 

0.  1C746 

0.C319C 

0. C2860 

0.  1C356 

0.03  161 

0.C2846 

0.C5576 

0.C311R 

C. 02817 

0.C96  39 

0.03065 

0.02777 

0.05383 

C.C30Q4 

0.02727 

0.C92  18 

0.C2934 

0.02666 

0.C9125 

0. 02852 

0.02596 

0.C8989 

0. C2749 

0.02513 

1.(6555 

0.  C2608 

0.02408 

0.07693 

0. C2386 

0.022  16 

0.C6298 

0. C 1988 

0.01891 

0.C4667 

0. C1345 

0.01293 

0. C3926 

0.CC4R5 

0.00462 

0.04597 

O.OOC72 

0. OCO68 

C . 05  2 0 1 

O.COOOR 

0.00009 

0.  17814 

2 Z 


TABLE  VI 


TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 435  km  Insteaci  of  385  km 


LCWi-:1  DENSITY  SATEL  1 F A tlF.CTOF  Y 155. SC  MHZ 


EANCE 

DELTA 

FANCE 

ELEVATION 

3 30*4 .9  2 

C . 40  36 

-0.  1533 

3117.18 

C.4024 

1.5267 

2929.60 

C. 3979 

3.3584 

2792.37 

C.  3898 

5.5244 

2555.76 

C.3781 

7. 4539 

2370.1  1 

C . 36  29 

9.7646 

2185.e5 

C . 3444 

1 2.  3664 

2003.53 

0.3231 

15.2659 

1823.94 

C . 2997 

18.5729 

1 6 • e . 15 

C.274R 

22.41C6 

1977.71 

C . 2492 

26.5466 

1314.66 

C. 22  36 

32.4202 

1162.59 

C.  1990 

39. 1371 

1027. 13 

C.  1765 

47.479C 

914.60 

C . 1577 

57.7878 

635.01 

C. 1442 

65.9579 

796.34 

0. 1380 

61.  16C5 

810.49 

0. 1401 

76.  1568 

666.39 

C. 1501 

63  .656  2 

966.49 

C. 1666 

52.32C7 

1C91.58 

0.1675 

43.0C62 

1236.  10 

C.2  11  3 

35.5141 

1393.91 

C . 2366 

29.4532 

156C.87 

C . 26  24 

24.4741 

1734.22 

C . 2878 

2C.3C5C 

1912. C9 

C.3121 

16.7460 

2C92.19 

C. 3345 

13.6523 

2276.60 

C.3543 

1C. 5163 

2461.68 

C.3711 

6.4670 

264  7.91 

0.3845 

6.2404 

2634.94 

C.3S41 

4.1553 

3022.45 

C.4001 

2.2563 

3210.22 

C.4026 

0.5238 

DELTA 

DELTA 

RELATIVE 

FI  . 

EL2 

EEFCR 

0.  C 2 9 4 3 

0.02553 

0.  1324  1 

0.03C93 

0.C2690 

C.  13042 

0.03184 

0.  C2786 

C.  12495 

0.  C 3 2 04 

0.02814 

0.  11544 

0.C3149 

0.02826 

0.  10257 

0.03020 

0.02756 

0.C6737 

0.C2830 

0.02628 

0.C7  1 15 

0.02591 

0.02448 

0.05519 

0. C2320 

0.  02225 

C.C4C67 

0.C2031 

0.01973 

G.  02828 

0.C1735 

0.C1703 

0.C  1858 

0.01442 

0.01425 

O.C1143 

0, 01155 

0.01147 

0.CC660 

0 , CC877 

0.00874 

0.C0384 

0.  00609 

0.C0608 

0.C0237 

0.00355 

0.00356 

0.CC  168 

0.00152 

0.00152 

0.0016  1 

0.  CO 240 

0.00239 

0.CC155 

0. 0C4B1 

0.C0430 

0.CC195 

0.CC744 

C. 00742 

0.C0290 

0.  01018 

0.01013 

0.  (C509 

0.01302 

0.01291 

0.CC9C1 

0 . C 1595 

0.01571 

0.C1506 

0.01893 

0. 01848 

0.02360 

0.C2188 

0.02111 

C. 03493 

0.02471 

0.02350 

0.04e66 

0. C2727 

0.C2552 

C.C6426 

0. 02943 

0.02705 

c . c ec  7 1 

0.03101 

0.02801 

0.C9668 

0.03  189 

0.C2836 

C.  11CP5 

0.C3201 

0.  028  11 

C.  122  05 

0.03138 

0. 02732 

0.  129  39 

0.C3010 

0.C2611 

0,  13258 
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TABLE  VII 


TEST  OF  ELEVATION  ERROR  FORMULA; 

H = 400  km  Instead  of  350  km 

NCFMO  DENSITY  FEE  NTEY  T E E J ECTOFY  155.50  H HZ 


DELTA 

DELTA 

DELTA 

RELATIVE 

RANGE 

RANGE 

ELEVATION 

EL. 

EL  2 

EFFOR 

3674.34 

1 .7C97 

0. 2219 

0.  14045 

0.11951 

0.  14912 

352C.74 

1.7045 

1.3635 

0.  14425 

0. 12289 

0. 146C3 

3364.98 

1.6530 

2.501 1 

C.  14660 

0.  12533 

0.14511 

32C7.C0 

1.6756 

3.6350 

0.  14758 

0.  12682 

0.  14062 

3046.75 

1.6528 

4.7654 

0.  14734 

0. 12747 

0.13436 

2 884.18 

1.6251 

5.8931 

0.  14611 

0. 12717 

C.  12824 

2719.23 

1.5933 

7.0183 

0.  144  14 

0. 12667 

C. 12117 

2551.6  5 

1.5574 

8.1415 

0.14167 

0. 12552 

0.11405 

2381.58 

1.5177 

9.2627 

C.  13896 

0. 12406 

0. 10724 

220  S . 56 

1.4734 

10.3820 

0.  13618 

0.  12242 

C . 10105 

2034.53 

1.4229 

11.4S9C 

0.13345 

0. 12069 

O.C9*56 

1656.64 

1.3639 

12.6130 

0.  13C80 

0. 11901 

0 . 090  1C 

1676.42 

1.2913 

13.7218 

0.  12804 

0. 11740 

C.C6305 

1493.23 

1.1540 

14. 8221 

0.  12440 

0. 11536 

O.C7265 

13C7.  19 

1.0472 

15.9C59 

0.11743 

0.11049 

0.05911 

1 118.28 

C . 6257 

16.9569 

0.  10304 

0.C9845 

0.04457 

926.54 

C . 5058 

17.9529 

0. C7384 

0.07104 

0.C37S2 

732.26 

C.  1873 

18.8427 

0.C3438 

0. C3270 

0.04876 

536.25 

C.C351 

19.5094 

0.00888 

0.00828 

0.(66  17 

34C.80 

C.0010 

19.3670 

0.00040 

0.00039 

0.04328 
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TABLE  VIII 


NOFH40 


TEST  OF  ELEVATION  ERROR  FORMULA; 
H = 400  km  Instead  of  350  km 


DENSITY  SATEL  TRAMCTCFY  155.50  MHZ 


RANGE 

DELTA 

FAKGE 

ELEVATION 

3306.21 

1 .6979 

-0.0839 

2118.47 

1.69C7 

1.6412 

293C.67 

1.6677 

3.4753 

274  3. 6 1 

1 .6280 

5.4406 

2556.56 

1.5718 

7.5663 

2371.25 

1.5CC3 

5.85C5 

2166.52 

1.4156 

12.4(37 

2004.53 

1.3203 

15.3533 

1624. 86 

1.2175 

18.6456 

1648.59 

1.1101 

22,  4765 

1476.46 

1.0013 

27.0C39 

1315.53 

C.6943 

32.4654 

1163.58 

C.7927 

39. 1729 

1C27.66 

C .7009 

47.5C58 

915. C7 

C.6245 

57.PC63 

835.44 

(.5704 

69.5666 

796.75 

(.5455 

81. 1651 

8 1C. 90 

C.  55  39 

76.2C40 

665.83 

C . 5943 

63.67(8 

966.98 

C.66C7 

52. 3434 

1092.14 

C .7459 

43.0395 

1236.73 

C.e437 

35.5547 

1394.62 

C . 9490 

29.5036 

1561.67 

1 .0578 

24.5351 

1735.10 

1.1668 

20.3767 

1913. C5 

1.2728 

16.8286 

2C94.22 

1.3724 

13.7453 

2277.71 

1.4628 

11.0207 

2462.85 

1.5409 

8.5769 

2649.13 

1 .6045 

6.3554 

2836.19 

1.6517 

4.3124 

3023.73 

1.6819 

2.4142 

3211.52 

1.6957 

0.6356 

DELTA 

DELTA 

RELATIVE 

EL. 

EL2 

EFROR 

0. 13875 

0.11757 

0 . 1 5 26  7 

0.  145  39 

0. 12349 

C.  15C67 

0. 14973 

0.12727 

0.14426 

0.  14824 

0.12B47 

0.  1 3340 

0. 14  339 

0. 12679 

c.  iieee 

C.  13611 

0.12222 

0.  10204 

0. 12561 

0.11502 

0.(8428 

0.  11326 

0.  10567 

C .(67C6 

C.C9990 

0.09476 

0.05147 

0.(8620 

0. C8291 

0.038  16 

0.C7267 

0.07068 

0.02750 

0.C5565 

0.05849 

0. C 1947 

0.C4728 

0.04663 

0.01371 

0.03560 

0.03524 

0.C1002 

0.C2458 

0.02439 

0.CC776 

0.01427 

0.01417 

0.00656 

0.00609 

0.00605 

0.CC614 

0.00962 

0.00956 

0.09631 

0. 01936 

0.01922 

0.CC701 

C. 030  10 

0.029e4 

C.CC87  1 

0.04152 

0.04103 

0.C1171 

0.C5363 

0.05275 

0.01645 

0.06644 

0.06489 

0.02337 

0.C7985 

0.C7723 

0.C2283 

0.C9362 

0.  C8941 

0.04500 

0.10731 

0.  10091 

0.(5970 

0.  12030 

0.  11 11  1 

0.C7641 

C.  13181 

0. 11939 

0.C9419 

0.  14C94 

0.  12520 

0.11173 

0.  14689 

0.  12815 

C.  12757 

0.14906 

0.12816 

0.  14026 

0.  14729 

0.12539 

0.14871 

0.  14  190 

0. 12027 

C.  15245 
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is  small,  then 


—A-  AR  « X(r)R 
so 

|[_^_  AR]  /R  | = X(r)R  (22) 

Since  N = 0 at  the  peak  of  the  profile,  the  height  h at  which  the 
derivative  (22)  changes  sign  can  be  taken  as  an  estimate  of  H. 
This  method  applies  primarily  to  post-mission  analysis,  but  also 
allows  correcting  a previously  used  H for  the  remainder  of  the 

re-entry  trajectory  in  real-time. 

It  is  also  possible  to  use  the  doppler  velocity  measured 

at  two  frequencies  to  generate  an  electron  density  profile  from 
which  H can  be  found.  To  see  this  we  next  consider  the  effect  of 
the  ionosphere  on  doppler  velocity. 

5 . DOPPLER  ERROR 

The  phase  of  the  radar  returned  signal,  relative  to  the 
local  oscillator,  is 

i . / nds  (23) 

C “'p 

where  oo  = 2TTf . Again  retaining  0(X)  in  n,  this  becomes 


5 = 2^  [R  - fx(s)ds] 

t J p 

• j 

The  rate  of  change  of  phase  is  § = $ > and  the  measured  doppler 

velocity  is 


VD  ~ Sw 


= Vt  ' tft 


f X(s) 


ds 


where  V,  = 


R^  is  the  true  doppler  velocity 


(24) 
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Combining  (2b)  and  (24) : 


VD  = Vt  - B/f2  (25) 

If  V is  measured  at  frequencies  f..  and  f„,  then  the  same 
D L & 

argument  as  used  for  range  gives  the  doppler  error  AV-^  as 

AV1  = vx  - vt  = g(V2  - Vx)  (26) 

where  g is  given  b/  (4) . 

It  should  be  noticed  that  if  the  velocity  V is  computed  by 
differentiating  the  range  in  (2a)  , the  result  is 

v_  = V.  + B/f2  (27> 

K L 

Thus  the  error  in  range  rate,  as  measured  by  time  delay,  is  the 
negative  of  the  error  in  coherent  doppler.  This  result  could  be 
employed  to  find  the  doppler  correction  using  a single  radar 
frequency : 


(28) 


The  doppler  correction  history  AV(t)  could  also  be  inte- 
grated backward  in  the  post-mission  analysis  of  a re-entry  to 
generate  the  range  correction  history  AR(t).  However,  this 
would  require  that  a hard  body  be  tracked  continuously  through 
the  ionosphere,  and  is  therefore  of  limited  usefulness.  Also, 
the  range  rate  VD  is  generally  not  as  precisely  determined  as 
V so  that  the  two-frequency  correction  (26)  is  preferable  to 
(28)  if  available. 


6.  ELECTRON  DENSITY  PROFILES 

Once  the  range  and  doppler  errors  have  been  found  over  an 
altitude  interval,  an  electron  density  profile  can  be  computed 
for  that  interval.  The  profile  will  of  course  be  of  electron 
density  along  the  trajectory,  rather  than  the  customary  vertical 
profile . 
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The  total  time  derivative  of  the  range  error  is  from  (1) 


Dt  ar  - |r  3r  + e + a aA  + Bt|  / X(s)ds  (29) 

where  E is  elevation,  A is  azimuth  , and  ==  . The  partial 

derivative  with  respect  to  time  results  from  secular  changes  in 
electron  concentration  along  the  ray  path.  Estimates  of  the 
maximum  rate  of  change  of  columnar  electron  density,  as  found 
from  satellite  measurements,  show  that  this  term  can  be  neglected 
for  the  applications  considered  here.  Also,  we  assume  that 
horizontal  gradient, >;  of  electron  density  are  small  compared 
to  vertical  gradients,  so  that  the  rate  of  change  of  R with 
azimuth  can  be  ignored.  Then  (29)  becomes 


DtAR  = X (R)R  + e 


X(s) ds 


(30) 


If  X(s)  is  assumed  to  be  a function  of  altitude  only,  we 
can  evaluate  the  rate  of  change  with  elevation  by  writing 


X(s , E + dE)  = X(s  + dp , E) 

where  dp  is  the  distance  which  must  be  added  to  s to  arrive  at 
the  same  altitude  as  that  resulting  from  an  increase  of  dE  in  E. 
From  Figure  3,  we  find  that 


dp  = s dE  cot  (E  + 8) 


= r s cos  E dE/(r  sin  E + s) 
o o 


Thus  the  derivative  becomes 
„R  -R 


J X(s)ds  = J 


[X(s  + dp ) -X(s) ]ds/dE 
R 


r 

= r cos  E / dX/ds  • s/(s  + r sin  E) 

° JL  ° 


ds 
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or 


R . r 1 2 1 

J X(s)  ds  = R • cot  E |X/(1  + B)  - J duX (u ) / ( 1 + Bu)  | 


(31) 


Using  equation  (16a)  for  6,  this  can  be  written  as 

R . I 

J X (s ) ds  = R/(l  + B)  ■ |X  cot  E - 6 | 


'o 


Combining  this  result  with  (30)  and  solving  for  X(R) : 
40.2  _ Df  AR  + 0^ 


(32) 


DtR  + q cot  E 


with 


q = R E/(l+B) 

In  (32),  the  quantity  R can  be  found  from  the  range  error. 
However,  it  is  preferable  to  use  the  error  in  coherent  doppler 
velocity  Z.Vp,  as  given  by  (26),  if  the  S/N  ratio  is  sufficient 
to  provide  accurate  doppler  velocities  at  both  frequencies.  If 
this  is  done,  it  must  be  remembered  that 

Dt  R = " VD  * 

In  (32) , all  quantities  except  the  elevation  error  6 are 
measured  directly  (6  could  also  be  measured  if  elevation  could  be 
measured  at  two  frequencies) . Using  (20)  for  6 implies  that  a 
peak  height  H for  the  ionosphere  must  be  chosen  in  order  to 
compute  the  profile.  In  practice,  this  can  be  handled  by  choosing 
a nominal  H,  computing  a profile  N,  and  using  this  profile  to 
determine  a revised  H.  If  desired,  this  process  can  be  con- 
tinued in  an  iterative  manner  to  obtain  successively  more  accu- 
rate profiles  and  H-values.  Alternatively,  if  the  look  angle  is 
small,  then  E is  small  and  (32)  can  be  approximated  by 


40.2 


N(h) 


D iR/D  R 


(33) 
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This  result  can  then  be  used  to  generate  a first  order  profile 
from  which  an  H can  be  determined  for  use  in  (32).  The  profile 
then  obtained,  if  not  very  different  from  the  first  order  re- 
sult (33)  , can  be  taken  as  correct. 

The  terms  in  (32)  were  estimated  for  relative  magnitude 
using  the  re-entry  data  from  Table  I.  It  was  found  that  at  the 
beginning  of  the  tabulated  trajectory  (E  < 10  ) , qt  was  com- 
parable to  D AR  and  q cot  E to  D^R  in  size.  Later  in  the  tra- 
jectory, at  altitudes  near  the  peak  height  H.  the  ratios  were 
q6/Dt.'R  ^ 0.06  and  q cot  E/D^R  « 0.18.  This  indicates  that  for  a 
typical  re-entry,  (33)  is  not  a good  approximation  at  long 
range,  but  is  fairly  good  near  the  peak  of  the  ionosphere.  It 
should  therefore  provide  a good  first  estimate  of  H. 

Once  a sufficiently  accurate  height  H and  density  profile 
is  found,  the  width  ah  can  also  be  determined.  If  AN  is  the 
full  3 db  width  in  terms  of  altitude,  then 

a,  = 0.6  AN 
h 

is  the  appropriate  value  to  use  in  (19).  However,  if  the  data 
are  insufficient  to  generate  a reliable  profile,  then  a fixed 
average  value  of  ah(e.g.  150  km)  may  be  used  in  (19),  since 
errors  in  a will  have  even  less  effect  on  6 than  errors  in  H. 

7.  VALIDATION  OF  ERROR  FORMULAS 

The  above  analysis  was  applied  to  a recent  mission  which 
contained  a long  track  on  a single  hard  target.  The  target  was 
tracked  from  above  700km  height  through  re-entry  by  ALTAIR  UHF 
and  VHF , and  from  above  600km  through  re-entry  by  ALCOR  C-band. 
The  two-frequency  formulas  were  applied  to  the  ALTAIR  data  to 
find  the  VHF  range  and  doppler  errors.  The  range  differences 
were  numerically  differentiated,  reversed  in  sign,  and  plotted 
as  a function  of  height,  along  with  the  doppler  differences. 
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The  results  are  shown  in  Figure  5.  The  agreement  between  -Dt'R 
and  V is  excellent,  and  the  two  plots  form  a well-defined 
profile.  Since  R is  almost  constant  over  the  interval  of  interest, 
the  profile  is  proportional  to  the  approximate  electron  density 
profile  given  by  (33) . Thus  Figure  5 was  used  to  obtain  a peak 
height  H = 340km.  This  H was  then  used  in  (20)  to  generate  the 
elevation  corrections  6. 

The  VHF  metric  data  (R,  AZ , El)  were  then  tested  against 
the  ALCOR  data,  both  with  and  without  the  ionospheric  cor- 
rections. To  do  so,  a least-squares  trajectory  was  fit  to  the 
ALCOR  data,  referred  to  ALTAIR  coordinates,  and  subtracted  from 
the  ALTAIR  data.  The  range  and  elevation  residuals  using  un- 
corrected VHF  data  are  shown  in  Figures  6 and  7.  The  residuals 

using  ALTAIR  data  with  ionospheric  correction  are  showi  in 

* 

Figures  H and  9. 

As  seen  in  the  figures,  the  range  error  without  correction 
is  as  large  as  1.7km  during  the  observed  part  of  the  trajectory. 
Applying  the  correction  (3)  reduces  the  maximum  error  to  -19m 
and  the  rms  error  to  -10m.  Elevation  error  without  correction 
reaches  0.17°.  Using  the  correction  (20)  reduces  the  maximum 
error  to  0.05°  and  the  rms  error  to  0.021°,  which  is  of  the  same 
order  as  the  expected  residual  elevation  error  of  the  ALTAIR 
system . 

These  results  can  be  taken  as  additional  validation  of  the 
error  formulas,  and  as  a demonstration  of  their  application  to 
metric  data. 


* In  each  case,  a constant  bias  of  0.14  has  been  subtracted 

from  the  ALTAIR  elevation  before  plotting.  This  bias  is  of  un- 
known origin  and  is  approximately  the  same  on  all  recent 
missions.  It  can  be  determined  by  comparison  with  ALCOR  ele- 
vation during  the  portion  of  a trajectory  lying  below  the 
ionosphere . 
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DIFFERENCE 


RANGE  RESIDUALS  ( m) 
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8.  FARADAY  ROTATION 


An  additional  error  in  measured  metric  parameters  is  caused 
by  the  earth's  magnetic  field.  Propagation  of  a linearly 
polarized  wave  in  the  ionosphere  through  the  earth's  field  re- 
sults in  a rotation  of  the  plane  of  polarization  of  the  wave. 

For  the  circular  polarization  transmitted  by  ALTAIR,  the  effect 
can  be  described  by  an  additional  correction  in  the  refractive 
index : 


n (s)  = 1 - X(s) 


1 + t,  (s) 


(34) 


with 


11  = 


where 


Wr 

— cos  $ 
GO 

eB 

me 


is  the  electron  gyrof requency  in  the  earth's  magnetic  field  B. 
and  ? is  the  angle  between  the  field  and  the  direction  of  propa- 
gation. In  (34)  T1  « 1 has  been  assumed,  along  with  quasi-longi- 
cudinal  propagation,  which  requires  that  $ - 90°  > 3°.  The 

upper  sign  in  (34)  applies  to  the  ordinary  ray  and  the  lower 


sign  to  the  extraordinary  ray. 


(3) 


The  range  error  resulting  from  (34)  is 


AR  = (l  + T|  | (£R)ion  (35) 

where 

TT  = <C°R  cos  i 

oo  I 

is  an  average  over  the  ray  path,  and  (AR)ion  is  the  range  error 

neglecting  the  earth's  field,  given  by  (1). 

To  estimate  the  size  of  this  correction  at  VHF , choose 

o 

u)g  8 • 10  /sec.  and  i * 60  for  a typical  re-entry,  so  that 
T\  s-  0.004.  For  propagation  northward  (along  the  magnetic  field, 
i.e.  $ < 90),  the  right  circularly  polarized  beam  transmitted  by 
ALTAIR  is  the  extraordinary  ray,  so  that  the  lower  sign  in  (35) 
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applies,  and  the  range  error  due  to  Faraday  rotation  is 

(AR)b  = 0.004  (AR)ion 

For  the  test  data  presented  above,  the  maximum  (AR^Qn  is 
1.7km,  resulting  in  (AR)B  « 7m.  The  largest  range  error  typi- 
cally seen  by  ALTAIR  is  about  3km  , corresponding  to  (AR)g  « 12m. 

Although  only  right  circular  polarization  is  transmitted, 
both  left  (principal)  and  right  (orthogonal)  circular  polari- 
zation is  received  by  ALTAIR.  The  right  circular  return  will 
show  no  net  magnetic  field  effect,  since  the  integrated  effect 
along  the  return  path  cancels  that  of  the  outward  path.  Thus 
at  a fixed  radar  frequency,  the  range  error  of  the  principal 
polarization  is  larger  than  that  of  the  orthogonal  polarization 
by  the  amount  (AR)  B . Therefore,  at  VHF , range  to  a target 
measured  at  left  circular  polarization  may  be  as  much  as  12m 
larger  than  the  range  measured  at  right  circular  for  a typical 
re-entry  trajectory.  This  range  difference  is  a significant 
percentage  of  the  30m  VHF  range  resolution.  Thus  it  should 
be  taken  into  account  when  detailed  comparisons  are  made  of  left 
and  right  circular  data,  as  in  the  study  of  sparse  chaff  clouds. 

Similarly,  the  earth's  field  causes  a difference  in 
doppler  velocity  measured  at  the  two  polarizations.  The  differ- 
ence is  ( A V ) „ = T (AV)  . , and  using  Tf  =0.004  and  (AV).  20m 'sec 

D i on 

for  the  ionospheric  error,  we  find  (AV)g  « 0.08m/sec.  Thus  the 
effect  of  Faraday  rotation  on  doppler  velocity  can  be  ignored 
for  most  applications.  Also,  the  effect  on  measured  elevation 
can  be  neglected,  since  it  is  generally  no  larger  than  0.4%  of 
the  ionospheric  correction  given  by  (20) . 

Further,  it  is  not  necessary  to  consider  the  effect  when 
computing  range  corrections  by  the  two-frequency  method.  To  see 
this,  one  recomputes  (3)  using  (34)  for  the  refractive  index. 
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The  result  is 


1 ‘l  TTp 


ARv  - (1  + g)  [Rv  - Ru] 

where  g is  given  by  (4)  and  p = f /f  « 0.375.  Thus  the  effect 
of  using  the  refractive  index  (34)  is  to  change  the  estimate  of 
the  VHF  range  error  by  a factor  (1  - 0.102  T|v)  . Again  assuming 
~ = 0.004  and  a maximum  value  of  3km  for  H Ry , the  additional 

range  correction  is  1.2m.  For  the  test  data  presented  above, 
the  correction  is  0.7m.  These  values  are  small  compared  to  the 
ALTAIR  system  range  residuals  and  can  therefore  be  neglected. 

In  summary,  the  errors  in  ALTAIR  metric  data  produced 
by  Faraday  rotation  are  negligible  when  two-frequency  measure- 
ments are  used  to  correct  for  ionospheric  refraction.  However, 
it  is  important  to  consider  the  effect  on  range  when  comparing 
data  taken  at  principal  polarization  with  that  taken  at  orthogo- 
nal polarization. 

9.  CONCLUSION 

The  ionospheric  corrections  described  here  are  currently 
being  employed  at  Lincoln  Laboratory  in  the  analysis  of  ALTAIR 
data.  As  demonstrated  in  the  test  data  presented  here,  these 
corrections  have  been  quite  successful  in  reducing  metric  residu- 
als, as  defined  by  reference  t:»  ALCOR  trajectory  data.  The 
corrections  have  also  proved  helpful  in  improving  impact  and 
pierce  point  predictions  based  on  exoatmospheric  tracks. 

The  corrections  can  also  be  applied  in  real  time,  although 
onlv  the  range  correction  (i)  has  been  used  systematically  thus 
far.  Adding  the  elevation  correction  (20)  would  result  in  im- 
proved trajectory  information  for  reacquisition  and  hand-off  to 
other  radars. 
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